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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace 
Remote Sensing program, AgRISTARS, is a program of research, develop- 
ment, evaluation, and application of aerospace remote sensing for 
agricultural resources. This program is a cooperative effort of the 
National Aeronautics and Space Administration, the U.S. Departments of 
Agriculture, Conmerce, and the Interior, and the U.S. Agency for Inter- 
national Development. AgRISTARS consists of eight individual pro- 
jects. 

The research reported herein was primarily sponsored by the 
Inventory Technology Developihv't (ITD) Project under the auspices of 
the Earth Resources Applications Division of the NASA/ Johnson Space 
Center. Dr. Jon Erickson is the NASA Manager of the ITD Project and 
Mr. Lewis Wade was the Technical Monitor for the reported effort. 

Included in this report is a summary of research sponsored by the 
Supporting Research (SR) Project under the auspices of the Earth Re- 
sources Research Division of the NASA/ Johnson Space Center and pre- 
viously documented within that project. Mr. Robert B. MacDonald was 
the NASA Manager of the SR Project and Dr. Glen Houston was the Technical 
Coordinator for the referenced effort. 

The analysis of corn and soybean profile characteristics was per- 
formed within the Envirormiental Research Institute of Michigan's 
Infrared and Optics Division, headed by Marvin R. Holter, Vice- 
President of ERIM, under the technical direction of Robert Horvath, 
Program Manager and Richard C. Cicone, Task Leader. 
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INTRODUCTION 


Central to the successful use of remotely sensed data for agri- 
cultural Inventories Is the ability of human analysts or computer 
algorithms to detect differences In the spectral characteristics of 
various cover classes. Experience with Landsat data In the Large Area 
Crop Inventory Experiment (LACIE) and other studies has demonstrated 
that the use of mul t1 temporal , spatially-registered data greatly en- 
hances the ability to distinguish between various crop spectral 
patterns [1]. In recent years Interest has been renewed In utilizing 
characterizations of the continuous patterns of crop spectral develop- 
ment over time, termed "profiles". In automatic crop Identification 
techniques [2,3,4]. These and other automated approaches offer sub- 
stantial gains In efficiency over manual techniques, and would therefore 
be of great value If their accuracies were similar to or better than 
those associated with human analysts. 

One substantial cause of error In both automatic and manual crop 
labeling techniques Is the deviation of a crop from Its expected 
spectral pattern due to cultural or environmental Influences (e.g., 
fertilization, moisture stress, changes In planting practices, etc.). 

This spectral deviation Is the result both of physiological changes In 
the plants themselves, resulting In changes In the spectral properties 
of the plant parts, and of changes In the canopy geometry. Including 
the orientation of plant parts, number and size of leaves or other plant 
parts, and amount of soil visible through the canopy. Particularly for 
computer algorithms, which lack the flexibility and adaptive capabilities 
of the human mind. It Is essential that major external Influences on 
crop spectral patterns be known In advance and taken Into account. If 
the algorithms are to perform adequately over broad regions or many 
growing seasons. 


gll 

Understanding the general patterns of crop spectral development, 
and the Influences of field conditions or those patterns. Is also essen- 
tial for spectrally-based assessment of crop condition. When combined 
with a means of estimating the time of occurrence of key stages of crop 
development, such a capability could contribute substantially to our 
ability to accurately estimate yields. 

This report presents the results of research aimed at characteriz- 
ing and understanding the spectral development patterns of corn and soy- 
beans, using field-collected reflectance data. Average profiles are 
described, as are the changes in those profiles brought about by changes 
In some major cultural and environmental factors. In addition, the 
association of profile features with stages of development of the two 
crops is discussed. Finally, the separability of corn and soybeans 
both In a general sense and in the context of particular field condi- 
tions Is considered. 

Development of the profile analysis technique described In Section 
2, and the literature review and Initial evaluation of cultural and 
environmental influences on Green Reflectance profile features which are 
Included in Section 5, were carried out under the auspices of the Sup- 
porting Research project, and were reported in that project [5]. Some 
of the results in Section 6 were previously presented in an ITD report 
[6]. The present report extends and summarizes the entire analysis. 
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DESCRIPTION OF DATA 

Evaluation of crop spectral characteristics can best be accom- 
plished with data collected at frequent Intervals over plots whose con- 
ditions are controlled or known. Such data have been collected for 
several years by ani at Purdue/LARS as part of a field research program 
carried out for NASA. For the analyses reported In this report, data 
were selected from experiments carried out in the 1978 through 1980 
growing seasons, which Included as experimental treatments nitrogen 
fertilization, planting date, and plant population for corn and variety, 
planting date, and row spacing for soybeans [7,8,9]. Reflectance mea- 
surements were made on clear days, resulting in gaps between successive 
observations of several days to weeks. All reflectance data were col- 
lected as or converted to Landsat-MSS inband reflectance values,* and 
multiple observations of a single plot on a single day were represented 
by their mean. 

A reflectance equivalent of the Tassel ed Cap transformation [10] 
was used to provide spectral variables that were physical ly-interpretable. 
This transformation, in Its Landsat-MSS form, captures 95* or more of 
the total data variability over agricultural areas In two variables. A 
rotation of the first two principle components In the reflectance data 
set was used to derive reflectance equivalents of the two variables [5], 
termed Green Reflectance and Bright Reflectance. These two variables 
contained 99* of the total variability in this data set. Green 


Note: After completion of these analyses, a preprocessing error 
(which occurred after receipt of the data from LARS) was discovered 
which invalidated the data collected using the Exotech 20C spectro- 
radiometer (one of two instrimients used). Of the experiments included 
in these analyses, only the com nitrogen experiments contained such 
data, and these experiments also included useable data collected 
using the Exotech 100 Landsat-band radianeter. The effects of the 
bad data on the analysis results appear to be negligible. 
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Reflectance serves as a green vegetation Indicator, while Bright Reflec* 
tance Is related to soil brightness or plot albedo. 

The total data set consisted of observations from 118 corn plots 
and 171 soybean plots. However, some plots were not suitable for all 
analyses, so the actual number used In any particular evaluation varied. 
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PROFILE ANALYSIS TECHNIQUE 

Although the field reflectance data set provided more frequent 
observation* and more detailed agronomic Information than wc.u1d be 
available using Landsat data* the problems of temporal gaps In spectral 
measurements and Inexplicable data variations were still apparent. As 
a result. It was necessary to devise a technique by which profile 
'alues could be Interpolated between actual measurements* and some 
smoothing of the measured values could be achieved. One approach to 
accomplishing the^;e objectifies Is use of a mathematical model to des- 
cribe the spectral development patterns [3]. Such models have been de- 
veloped for and successfully applied to spring small grains [11*12]. 
However* the plateau feature commonly observed In corn Greenness and 
Green Reflectance data (Figure 1) cannot be adequately described by 
previously developed profile models [5], so some other approach •• 
either a new model or a more general technique - was required. 

Two techniques were selected for profile characterization. For 
analysis of soybeans, which exhibit a less complex profile shape In 
Green Reflectance, and for some analyses of corn* a cubic smoothing 
spline [13] was used. This method provided the desired degree of 
smoothing, and captured the plateau feature given that sufficient data 
points were present In that portion of the profile. However, when 
significant data gaps existed In critical periods, the cubic smoothing 
spline result for corn was less acceptable. Acco»*d1ngly , a profile 
model form, developed at ERIM specifically for this purpose [5], was 
used on the corn data when the plateau feature was of prime Interest. 

For Bright Reflectance data* early season variations caused by 
soil mo1stur“e differences and other unknown factors rendered automatic 
curve fitting Impractical. Figure 2 Illustrates the soil moisture 
effect for bare soil plots. Therefore* B.'lght Reflectance profiles 
were derived for each plot manually, based on the available data points 


5 





Bright Reflectance 


2pi 


ORIGINAL PAGZ IS 
OF POOR QUALITY 


Corn Experiments - Bare Soil Plots 
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and soil moisture information. Since planting dates were known, a "days 
:,ince planting" time axis was used in describing the comparing both 
Green Reflectance and Bright Reflectance profiles. 

In order to facilitate comparison of profiles, a set of features 
was defined which described their major characteristics. These features 
(illustr .ed in Figure 3) included the peak profile value (PMAX), the 
time at which the peak occurred (PT), the times of occurrence of one- 
half the peak value (HPl and HP2), and the time Intervals between HP1 
and PT (SPANl). PT and HP2 (SPAN2), and HPl and HP2 (SPAN3). In Green 
Reflectance, these features are related to the maximum amount of green 
vegetation (PMAX), the rate of vegetative development (PT, HPl, SPANl), 
the rate of senescence (SPAN2), and the overall development rate (HP2 
and SPAN3). For corn, two additional features id«re used to describe 
the duration and slope of the plateau (Figure 4). 

Features of Green Reflectance profiles were derived explicitly, 
and analyzed quantitatively by analyses of variance. Fo*' Bright Re- 
flectance, the features were derived implicitly by visual analysis of 
the profiles, and compared qualitatively. 
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FIGURE 4. ADDITIONAL CORN PROFILE FEATURES 
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NORMAL PROFILES 



Figures 5 through 8 illustrate the average profiles of com and 
soybeans, using the corn profile model results for 57 corn plots and 
the cubic smoothing spline results for 167 soybean plots. The dashed 
lines around the Green Reflectance profiles represent one standard de- 
viation about the mean. The asymmetry of the soybean variability about 
the mean, most noticeable in the declining phase of the profile, is 
largely the result of differences in timing of leaf senescence between 
maturity classes. 

The Bright Reflectance profiles, illustrated with a dark and light 
soil background, are delimited by dashed lines representing the range 
of mean profiles (by treatment) used to derive the overall average pro- 
files. It should be noted that for corn, the very bright soil com- 
pletely oDscures the peak in the Bright Reflectance profile seen with 
darker soils. For both corn and soybeans, soil effects were undetect- 
able in the latter portions of the profiles. 

Figures 9 and 10 present the spectral trajectories of the two crops 
in Green Reflectance and Bright Reflectance, using the dark soil Bright 
Reflectance profiles. The plateau in corn Green Reflectance, when com- 
bined with the Bright Reflectance profile, is expressed as a movement 
away from and then back to the "Green Arm" of the Tassel ed Cap. Figures 
11 and 12 present the same trajectories using the bright soil profiles, 
and illustrare the substantial effect of soil reflectance on crop trajec- 
tories in the early portion of the growing season. 
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CULTURAL AND ENVIRONMENTAL EFFECTS ON PROFILE FEATURES 

This section sumnarlzes the results of both a literature review of 
the effects of environmentel factors on corn and soybean character- 
istics and the profile Teci‘»res analyses. Reference [5] contains a more 
detailed description of many of the results. For Green Reflectance pro- 
files, results reported are those found to be significant to the 0.9 
level of confidence. Figures Illustrating the effects of axperlmental 
treatments are Intended to be descriptive rether than quantitative. 

5.1 CORN EFFECTS 

Nitrogen Fertilization . The availability of nitrogen, which Is 
required for synthesis of chlorophyll. Influences the vegetative de- 
velopment of corn. Abundant Nitrogen results In more and larger leaves, 
longer vegetative stages, and Increased longevity of green leaf area. 
These effects were expressed In a later and higher peak Green Reflectance 
value and a longer and flatter plateau (Figure 13a). While the quality 
of Bright Reflectance data for these experiments was low, some indi- 
cation of a higher peak Bright Reflectance value was observed (Figure 
13b). 


Planting Date . Later-planted corn experiences higher temperatures 
at any given stage of development than does earlier -pi anted corn. As a 
result, the rates of plant emergence arnl development are Increased, as 
are the rates of leaf emergence and leaf area development. Very early 
planting exposes the plants to lower temperatures, and would thus be 
expected to delay emergence and retard early growth. 

Both late and very early planting reduced the peak Green Reflec- 
tance profile value (Figure 14a), probably an Indication of the less 
conducive growing conditions encountered. Both also reduced the over- 
all development time (HP2). Late planting. In addition, hastened the 
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FIGURE 13. NITROGEN FERTILIZATION EFFECTS ON CORN PROFILES 
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rate of green-up, the time of peak, and the rate of green decline. 

Late planting similarly hastened the time of occurrence and reduced the 
value of the Bright Reflectance profile peak (Figure 14b). 

Plant Population . Increases In corn plant population density 
cause reduced rates of leaf area production, faster early height In- 
crease, greater ma>1mum leaf area Index, and faster decline In LAI after 
peak. These effects were expressed In an earlier and higher peak Green 
Reflectance value, a faster rate of green-up, and a steeper plateau 
slope (Figure 15a). In Bright Reflectance, a higher peak value was ob- 
served, and the elimination of soil effects occurred earlier. Indicat- 
ing earlier canopy closure (Figure 15b). 

5.2 SOYBEAN EFFECTS 

Variet y. Soybean varieties exhibit differences in days to maturity 
(maturity class), plant height, leaf size, number, and orientation, 
rate of accumulation and maximum leaf area, ability to achieve full 
closure, response to row spacing or planting delays, etc. Not surpris- 
ingly, all the profile features were significantly affected by varietal 
differences at some row spacing (Figure 16). 

Planting Date . Both early and late planting tend to cause reduc- 
tions in the final height of soybean plants, and to reduce the rate of 
canopy closure. The higher temperatures associated with later planting 
hasten emergence and early growth, and reduce the duration of the 
vegetative phase. 

The peak Green Reflectance value was reduced with both early and 
late planting, and was substantially earlier for late planting (Figure 
17a). In addition, late-planted soybeans had a faster rate of green-up 
and a shorter overall development time. Similarly, late-planted soy- 
beans exhibited an earlier peak Bright Reflectance value and a more 
rapid increase in Bright Reflectance, as well as a reduced overall pro- 
file span (Figure 17b) . 
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Row Spacing . Soybean varieties differ In their response to row 
spacing. In general, however, wider rows tend to cause a reduction In 
the rate of leaf area accumulation and a delay In achievement of full 
canopy closure. These effects were expressed In lower and later peak 
Green Reflectance values and a slower rate of green-up (Figure 18). 
Also observed was a more rapid green decline after peak, probably an 
Indication of the lower leaf area density, which allowed brown lower 
leaves or soil to show through the canopy sooner after the peak. 
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ASSOCIATION OF SPECTiyVL AND DEVELOPMENTAL EVENTS 

Another aspect of understanding crop profiles Is knowing the 
stages of development associated with some of the key profile features. 
Such knowledge allows development of crop Identification techniques 
which utilize the most fundamental and therefore most stable differences 
between crops, and could aid the assessmei.w of crop conditions or pre- 
diction of yield by providing a means of pinpointing certain key de- 
velopmental events from spectral data. 

In order to carry out this analysis, development stage data col- 
lected for the experimental plots were smoothed and Interpolated by 
polynomial regression. The time of occurrence of each stage, or the 
stage at any particular time, could then be easily determined. Figure 
19 Illustrates the result of combining the development stage data with 
the Green Reflectance data for a typical corn plot. 

6.1 CORN RESULTS 

Stages of development defined by Hanway [14] were used for corn. 
Peak Green Reflectance was found to occur at stages 2.5 to 3.0, which 
correspond to 10 to 12 leaves fully emerged. These stages occur about 
two weeks prior to tassel emergence, and three weeks before the stage 
normally associated with peak LAI (stage 5 - silking). Percent cover 
data associated with the experimental plots suggest that the Green 
Reflectance peak occurred before maximum canopy closure. Although 
this result seems to contradict the normally expected correlation 
between vegetation indices and LAI or percent cover, a plausible, if 
hypothetical, explanation can be given. 

At stage 2.5 to 3.0, all or nearly all the green leaf area is de- 
veloped, but much of it is still furled into a pseudostem (Figure 20a). 
Because the stem itself is only about half as tall as the total corn 
plant at this point, there is a dense and fairly shallow layer of pure 
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green leaf matter at the top of the canopy » which contains all the 
green leaf matter that the plant will have. Development after this 
point, Including continued stem elongation, changes In leaf angular 
orientation, and emergence of tassels, all serve to reduce Green Re- 
flectance by Increasing shadowing or changing the mix and distribution 
of plant parts In the c&nopy (Figure 20b). 

The plateau in corn Green Reflectance ends around stage 8, the 
early dent stage. This stage falls about one month after Initiation 
of rapid dry matter accunxjlatlon In the kernels, which would correspond 
to an accelerated rate of senescence In the vegetative parts. Since 
one would expect a time lag between onset of senescence, which proceeds 
from the bottom of the plant to the top, and any noticeable effect on 
the Green Reflectance of the canopy, the observed delay is not surprising. 

6.2 SOYBEAN RESULTS 

The stages of soybean develo|»nent defined by Fehr and Caviness 
[15] were used. Since many varieties of soybeans are Indeterminate, the 
Fehr/Caviness system provides separate vegetative and reproductive stage 
progressions. 

Peak Green Reflectance occurred between vegetative stages 12 and 
21 (12 to 21 nodes with fully emerged leaves), or reproductive stages 
3.5 to 6 (beginning pod to full seed). The extremes of these stage 
ranges normally occur about 30 days apart, suggesting that there Is 
little If any correlation between peak Green Reflectance and any par- 
ticular soybean stage of development. 

In some of the plots, peak Green Reflectance did coincide with the 
maximum vegetative stage reached, but In many others, the Green Reflec- 
tance peak occurred at a vegetative stage well before the maximum. In 
most of these cases, nowever, lodging was reported at the time of 
maximum vegetative development, and at all observation times after the 
peak. While the severity of lodging was not recorded, one can speculate 
that It was enough to reduce the Green Reflectance of the canopy by 
exposing more stems and changing the overall geometry. 
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It would appear, then, that maxlimim Green Reflectance In soybeans 
occurs at the end of vegetative develc^xnent except when lodging occurs. 
When lodging was a factor In these data, the end of vegetative develop- 
ment sometimes occurred at a point well down the declining side of the 
Green Reflectance profile^ Thus It must be concluded that Green Reflec- 
tance profile features, by themselves, can give little or no reliable 
Information as to the stages of soybean vegetative or reproductive 
development. 

If by some means It could be determined that a particular field, 
or a region In general, was planted with a determinate variety and/or 
that no lodging was present, then estimates of stage Information might 
be extractable from the profile features. However, since In most regions 
one could expect to find both determinate and Indeterminate varieties 
of several maturity classes and a range of susceptibilities to lodging, 
the likelihood of obtaining such Information for any given field se«ns 
slight. 
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SEPARABILITY OF CORN AND SOYBEANS 


Based on the profile-derived features for the entire set of ex- 
perimental plots, an analysis of the separability of corn and soybeans 
profiles was carried out. Both the corn model and cubic smoothing 
spline were used to characterize corn Green Reflectance profiles, the 
model because It more accurately describes the actual spectral develop- 
ment pattern, and the spline because It more closely resembles the kind 
of approach that might be used In an operational setting where the crop 
type of a sample w'^s not known. Each feature was hlstogrammed, and the 
histograms were compared to determine separability. Evaluation of 
Bright Reflectance profile separability was carried out qualitatively. 

7.1 OVERALL RESULT 

Using the corn model, the peak Green Reflectance profile value 
(PMAX), the time of occurrence of that peak (PT), the rate of vegetative 
development after emergence (SPANl), and the rate of Green Reflectance 
decline after the peak (SPAN2) all provided substantial separability. 

Corn tended to reach a lower peak value earlier, had a more rapid re- 
lative green-up rate, and declined In Green Reflectance much irore slowly 
than soybeans. ”e separability related to green decline was the result 
of the corn Green Reflectance plateau. 

When the cubic smoothing spline was used for both crops, substantial 
separability was still found In the height of the Green Reflectance pro- 
file peak (PMAX), with reasonably good separation also In the rate of 
Green Reflectance decline (SPAN2). Most notably, In this data set lOOX 
separability was achieved using the peak Green Reflectance profile value 
and the rate of Green Reflectance profile decline (Figure 21). While 
the peak value, or something similar to It, Is a feature used In several 
current corn/soybean discrimination techniques [16], the plateau feature 
Is little used at this time. However, apparent rates of green-up and 
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decline, as derived by simpler profile models, are used In current 
techniques [17], and the separability found In these features Is prob> 
ably an Indirect result of the corn plateau. 

Early season features such as the rate of early vegetative develop- 
ment (HPl) provided little or no separability, particularly with the 
spline technique. In addition, comparison of the slopes of the ascend- 
ing portions of the average corn and soybean profiles revealed little 
or no difference between the two crops In this respect (Figure 22). 
Planting date differences. Ignored In this analysis, may provide a 
greater potential for separation, depending on local crop calendars, 
but there Is little indication of purely spectral separability In the 
early season (I.e., before peak Green Reflectance). 

In Bright Reflectance the only obvious source of separability was, 
again, the height of the profile peak. As can be seen by comparing 
Figures 9 and 10, this is an expression of the same phenomenon expressed 
in the Green Reflectance profile peak height - Soybeans moves farther 
up the Green Arm" of the Tassel ed Cap than does corn. This feature 
too has been used In crop Identification techniques [18]. 

7.2 EFFECTS OF CHANGES IN FIELD CONDITIONS 

Comparison of the results In Section 5 and 7.1 reveal that many of 
the environmental and cultural factors considered affect precisely those 
features most Important In corn and soybean discrimination. Both 
nitrogen fertilization and increased planting density tend to raise the 
peak Green Reflectance value of corn, while early or late planting and 
wider row spacing tend to lower the peak soybean Green Reflectance 
value. Similarly, nitrogen deficiency and late planting tend to shorten 
or soften the plateau effect in corn. Under particular sets of con- 
ditions, then, one should expect a degradation In the separability of 
corn and soybeans. 
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CONCLUSIONS 

Comparison of the curves fit to each plot in the set of corn end 
soybeans transformed field ref1ectaii\.e data indicates that there are 
indeed characteristic profile shape^; for these two crops, shapes which 
show some variation but which nonetneless retain characteristic attri- 
butes. At least some of the variation is caused by changes in field 
conditions and cultural practices, and it has been shown that variations 
in these factors within a r^^ge that could be reasonably expected in an 
operational setting caure r grificant changes in profile features. 
Furthermore, while corn and soybeans are distinguishable based on their 
profile features, t>e effects of changes in conditions or practices are 
such that separability could, under particular sets of circumstances, 
be substantially degraded. 

The evaluation of development stage association with profile fea- 
tures shows a clear relationship between corn Green Reflectance features 
and stages of development, and further shows that the peak in corn Green 
Reflectance occurs earlier than would have been expected, before either 
maximum LAI or maximum canopy closure. Conversely, soybean Green Reflec- 
tance features cannot be associated with any particular stage of de- 
velopment. In the absence of lodging, a strong correlation can be seen 
between peak soybean Green Reflectance and maximum vegetative develop- 
ment. 


The profiles described in this report, and the changes in those 
profiles attributable to field conditions and cultural practices, pro- 
vide a foundation for development of automatic crop identification 
techniques using Landsat data. If, using meteorological or other avail- 
able information, field conditions could be inferred for a particular 
region of interest, expectations regarding corn and soybean profile 
shapes could be modified based on the presented results, an approach 
which could yield a substantial increase in labeling accuracy. A 
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technique adaptable in such a manner, and therefore applicable over 
wide regions and many years, would be of considerable value. 

With regard to assessment of field conditions the results des- 
cribed in this report suggest that spectral data alone may not provide 
adequate information for determining the presence or absence of a par- 
ticular physiological stress. For example, late planting and nitrogen 
deficiency, which primarily influence the physiological development of 
the individual plants, reduce the peak Green Reflectance value, but a 
similar result is seen when population density is reduced. Reduced 
plant population will affect individual plant growth to some degree, 
but its primary influence on the Green Reflectance profile is probably 
due to changes in the proportion of soil background as compared to vege- 
tation or the density of the vegetation viewed by the sensor. 

Although a thorough analysis of all profile features could allow 
discrimination between these various factors, it is most likely that 
the acquisition intervals provided by Landsat, and the potential for 
mixing several such factors in any given field, would only allow use 
of more gross characteristics such as the peak value. 

The likely impact of this result on spectral contributions to 
yield estimation is uncertain, however, since the final goal is not 
accurate identification of the yield-affecting condition but rather 
accurate estimation of its effect. In the previous example, reduced 
population dinsity may not be classified as a "stress", but it most 
probably will reduce the yield obtained from the particular field. In 
addition, the strong association observed between the Green Reflectance 
profile peak in corn and a particular stage of development suggests that 
a "critical time interval" for yield-affecting stresses could be esti- 
mated for each field. With this information, the likely impact on yield 
of stresses detected or forecast by other means could be more accurately 
predicted, allowing more accurate prediction of regional production. 
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